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Abstract. We discuss the time evolution of dust-to-gas 
mass ratio in the context of multi-phase model of inter- 
stellar medium. Phase transition of interstellar gas is con- 
sidered to occur on a timescale of ~ 10 7-8 yr, according 
to a nonlinear open system model of interstellar medium. 
Since the phase transition changes the dust formation and 
destruction rates, the dust-to-gas ratio also fluctuates on 
the same timescale. This explains the scatter of the dust- 
to-gas ratios of spiral galaxies quantitatively, though we 
should note the large observational uncertainty. 
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1. Introduction 

Recent chemical evolution models of galaxies including the 
dust content are successful in explaining the dust amount 
of nea rby g alaxies (W ang 



1991: Liscnfcld & Ferrara 



1998 



Dwek [L998| Hirashita |l999a| here after H99; see also Tak- 
agi, Arimoto, & Vansevicius 1999 ). Supernovae (SNe) are 
the dom inant source of the formation of dust grains (Dwek 
fc Sc alo 1980| ) , and SN s hocks destroy grains (Jones et al. 
199J; Borkowski & Dwck |l995|) . Thus, in those models are 



dust content connected with star formation histories. 

In our previous work, H99, the dust-to-gas ratio was 
expressed as a f unction of metallicity (see also Lisenfeld 
& Ferrara 1998), which is also related to star formation 
histories. It confirmed the suggestion proposed by Dwek 
( |1998| ) that the accretion process onto preexisting dust 
grains is efficient in spiral galaxies. However, since the 
accretion is effective in cold clouds, the global efficiency 
of the accretion depends on the fraction of the gas in cold 
phase (Seab |l987t McKee |1989| ; Draine h.990[ ). Thus, the 
efficiency varies on a timescale of ~ 10 7-8 yr by the phase 
transition of the ISM (Ikeuchi |1988| ; McKee |198S| ). 

In this Letter, we combine the framework of H99 with a 
theoretical work on multiphase ISM by Ikeuchi & Tomita 
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(1983), whose limit-cycle model of the ISM phase transi- 
tion is applied to the result in Tomita, Tomita, & Saito 



( p9l| by Kamaya & Takeuchi ( p9?| ; hereafter KT97). 



In the limit-cycle model, mass fraction of each phase os- 
cillates continuously because of mass exchange among the 
components of the phases. The timescale of the phase tran- 
sition in the model is determined by a few parameters 
intrinsic to a spiral galaxy; sweeping rate of SN shocks, 
evaporation rate of the cold gas, and the cooling rate of 
the gas heated by SN shocks. Actually, a static solution 
as well as the limit-cycle solution for the filling factors of 
the three components is possible depending on the param- 
eters. However, the oscillatory behaviour (i.e., the limit- 
cycle model) of the filling factors is supported observa- 
tionally. Indeed, the observed scatter of the far-infrared- 
to-optical flux ratios of spiral galaxies (Tomita, Tomita, & 
Saito 1996) is interpreted through the limit-cycle model in 
KT97. KT97 suggested that the fraction of the gas mass 
(i.e., the mass filling factor) in the cold phase changes in 
the range of 0.1 to 0.7 (or more) on the timescale of 10 7-8 
yr (see also Korchagin, Ryabtsev, & Vorobyov 1994). 

This Letter is organized as follows. In §2, we investigate 
the variation of the dust-to-gas ratio due to the phase 
transition of ISM in spiral galaxies. Finally, we discuss 
the result in S3. 



2. Grain growth in multiphase interstellar 
medium 

The ISM in a spiral galaxy is composed of multiphase 
gas. McKee & Ostriker ( 1977 ) constructed the model of 
the ISM with three components in a pressure equilibrium: 
the cold phase (T ~ 10 2 K and n ~ 10 cm -3 ), the warm 



phase (T ~ 10 4 K and n ~ 10 1 cm 3 ), and the hot phase 
(T ~ 10 6 K and n ~ 10~ 3 cm" 3 ). Since the mass of the 
hot component is negligible in a galactic disc compared 
with the others, we only consider the warm and the cold 

gas. 

Dwek ( [1998 ) and H99 showed that the accretion of 
heavy elements onto preexisting dust grains is the domi- 
nant process for the growth of the dust content in spiral 
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galaxies. Thus, we concentrate on the effect of the phase 
transition on the accretion. 

The timescale of the grain growth through the accre- 
tion of heavy element, r grow , can be estimated by the du- 
ration of the collisions between heavy-elements atom and 
grains. According to Draine (1990), r grow ~ 5 x 10 7 yr in 
cold gas. Here, we should note that the accretion process 
is more effective in denser environments. The efficiency 
of the accretion is proportional to the square of the gas 
density if metallicity and dust-to-gas ratio are the same, 
since the densities of both metal and dust contribute to 
the efficiency. Therefore, among the three components of 
the ISM, we only consider the accretion process in the cold 
gas, the densest component of the ISM. 

According to H99, the increase rate of dust mass by the 
accretion process in a galaxy, [dM^/ dt] acc , is expressed as 



dM d 
dt 



VM g (l - f) 



(1) 



where V is the dust-to-gas mass ratio, M g is the total mass 
of ISM in the galaxy (i.e., = DM g ), f is the fraction of 
the metal in dust phase, and r acc is the accretion timescale 
of heavy elements onto dust grains (see Eq. 3 in H99). We 
note that the newly introduced parameter r acc is different 
from r grow , since r acc is the accretion timescale averaged 
over all the ISM phases. As commented above, the dust 
in the cold gas dominantly contributes to the accretion 
process. Thus, [dAId/dt] acc can also be expressed in the 
following way: 



dMd 
dt 



VX cold M g (l - /) 



(2) 



' grow 



where X co \d represents the mass fraction of the cold phase 
to the total mass of ISM. Here, we assume that the values 



of V and / are constant for each phase. McKee ( 1989 ) 
showed that the mixing of phases makes the difference 
in the V values among phases negligible. We also expect 
that / is treated as con stant for all phases because of 
the mixing (Tenorio-Tagle 1996). Combining equations ([j]) 
and (||), we finally obtain 

(3) 



A co ld 



According to Ikeuchi (198S) and KT97, A co id can vary 



with the range of 0.1 ^ A co id ^ 0.7 in 10 7 8 yr. Therefore, 
from equation (||), we see that r acc varies in the range of 



1 At < t 

-L-^'grow rO ' accr 



^ 1 0r grow on that timescale. 



3. Discussions 

We have shown in the previous section that the param- 
eter r acc , the typical timescale of accretion of heavy ele- 
ments onto dust grains, changes on a timescale of 10 7-8 
yr through phase transition of ISM. The range of r acc is 



estimated as 1.4T grow < t 



a , rr < 10T grow , which is typi- 



the parameter /3 acc (proportional to the efficiency of the 
accretion of heavy elements onto preexisting dust grains), 
defined in H99, changes by nearly an order of magnitude. 
Moreover, the timescale of the variation is much shorter 
than the typical timescale of the gas consumption in a 
galactic disc (<; I Gyr; Kennicutt, Tamblyn, & Congdon 
1994). Thus, the dust-to-gas ratio in a spiral galaxy ex- 
periences a short-term (~ 10 78 yr) variation with the 
amplitude of an order of magnitude. 

The short-term variation can be tested by examining 
nearby spiral galaxies. The dust-to-gas ratios of the spi- 
ral galaxies shows scatter around their mean values even 
if the meta llicity is almost the same (Issa, MacLaren, & 
Wolfendale 199C; see also H99). According to Figure I in 
H99, the theoretical lines almost reproduce the observed 
values. However, the dust-to-gas ratios of the Galaxy and 
M3I differ by several times. Both the galaxies lie in a 
range of 5 ^ /3 acc ^ 20. This means that we can explain 
the dust-to-gas ratios of these galaxies if /3 acc changes by 
more than 4 times. Indeed, the discussion in §2 demon- 
strated that /3 acc can change by more than 7 times on 
~ 10 7-8 yr because of the phase change of ISM. Thus, it 
is possible to explain the scatter of the dust-to-gas ratios 
of spiral galaxies by considering the phase transition. 

As for dwarf galaxies, we need another way to approach 
them, since the heavy element accretion in dwarf galax- 
ies is much less efficient than spiral galaxies due to their 



small metallicity (Hirashita 1999b). Because of their shal 



low gravitational potential, the mass outflow (e.g., Mac 
Low & Ferrara 1999] ) can be responsible for the dust-to- 
gas ra tio spread, as emphasized by Lisenfeld & Ferrara 
( |1998|) . 

We only have considered the dust formation process. 
However, we should also consider dust destruction. The 
dominant dust destruction occurs in the warm and hot 



cally 7 x 10 7 yr ^ T accr ^ 5 x 10 s yr. This means that 



phases in which SN shock waves propagate (Seab 1987 ). 
This means that the destruction efficiency is expected 
to show anticorrelation with A co id. If a galaxy is in a 
higher- X co id state, the dust destruction is more inefficient 
whereas the dust growth is faster. Thus, the variation of 
dust-to-gas ratio may become larger if we take into ac- 
count the dust destruction. 

Finally, we should note that it is still probable that 
the scatter is caused by observational uncertainty, since 
the dust-to-gas ratio is not a direct observable. However, 
from the discussion in §2, we can still propose that the 
dust-to-gas ratio varies on a timescale of ~ 10 7-8 yr by 
nearly an order of magnitude. 
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